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Abstract 

A detailed study of the Counter-Rotating Model (CRM) for generic electros- 
tatic (magnetostatic) axially symmetric thin disks without radial pressure is pre- 
sented. We find a general constraint over the counter-rotating tangential velocities 
needed to cast the surface energy-momentum tensor of the disk as the superposi- 
tion of two counter-rotating charged dust fluids. We then show that this constraint 
is satisfied if we take the two counter-rotating streams as circulating along elec- 
trogeodesics with equal and opposite tangential velocities. We also find explicit 
expressions for the energy densities, electrostatic (magnetostatic) charge densities 
and velocities of the counter-rotating fluids. Three speciflc examples are conside- 
red where we obtain some CRM well behaved based in simple solutions to the 
Einstein-Maxwell equations. The considered solutions are Reissner-Nordstrom in 
the electrostatic case, its magnetostatic counterpart and two solutions obtained 
from Taub-NUT and Kerr solutions. 

PACS numbers: 04.20.-q, 04.20.Jb, 04.40.-b 



1 Introduction 

Stationary or static axially symmetric exact solutions of Einstein equations describing 
relativistic thin disks are of great astrophysical importance since can be used as models 
of certain stars, galaxies, accretion disk and universes. Theses were first studied by 
Bonnor and Sackfield [|I|, obtaining pressureless static disks, and then by Morgan and 
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Morgan, obtaining static disks with and without radial pressure 0, 0]. In connection 
with gravitational collapse, disks were first studied by Chamorro, Gregory and Stewart 
Disks with radial tension have been also recently studied 0. Several classes of exact 
solutions of the Einstein field equations corresponding to static and stationary thin disks 
have been obtained by different authors ^-[1^, with or without radial pressure. 

In the case of static disks without radial pressure, there are two common interpreta- 
tions. The stability of these models can be explained by either assuming the existence 
of hoop stresses or that the particles on the disk plane move under the action of their 
own gravitational field in such a way that as many particles move clockwise as counter- 
clockwise. This last interpretation, the "Counter- Rotating Model" (CRM), is frequently 
made since it can be invoked to mimic true rotational effects. Even though this interpre- 
tation can be seen as a device, there are observational evidence of disks made of streams 



of rotating and counter- rotating matter |17, 18 



Disklike sources in presence of electrogmatic fields, specially magnetic fields, are also 
of astrophysical interest mainly in the study of neutron stars, white dwarfs and galaxy 
formation. In the context of general relativity models of disks for Kerr-Newman metrics 
static axisymmetric spacetimes with magnetic fields [|19] and corformastationary 



metrics [^, have been considered recently. Following the Ref. fSQ] the resultating 
disks can also be interpreted either as rings with internal pressure and currents or as 
two counter-rotating streams of freely moving charged particles, i.e. which move along 
electrogeodesics (solution to the geodesic equation in the presence of a Lorentz force). 

The aim of this paper is to perform a detailed study of the CRM for generic electros- 
tatic (magnetostatic) axially symmetric thin disks without radial pressure. The paper 
is organized as follows. In Sec. 2 we present a summary of the procedure to obtain thin 
disks models as rings with a purely azimuthal pressure and currents using the well-known 
"displace, cut and reflect" method extended to solutions of Einstein- Maxwell equations. 
In particular, we obtain expressions for the surface energy-momentum tensor and the 
electrostatic (magnetostatic) current density of the disk. Next, in Sec. 3, the disks 
are interpreted in terms of the CRM. We find a general constraint over the counter- 
rotating tangential velocities needed to cast the surface energy-momentum tensor of 
the disk as the superposition of two counter-rotating charged dust fluids. We then 
show that this constraint is satisfied if we take the two counter-rotating streams as 
circulating along electrogeodesics with equal and opposite tangential velocities. We also 
find explicit expressions for the energy densities, electrostatic (magnetostatic) current 
densities and velocities of the counter-rotating fluids. In following section. Sees. 4, 
three specific examples are considered based in simple solutions to the Einstein-Maxwell 
equations. The considered solutions are Reissner-Nordstrom in the electrostatic case, 
its magnetostatic counterpart, and two solutions generated from Taub-NUT and Kerr 
solutions. In particular, we study the tangential velocities, mass densities and electros- 
tatic (magnetostatic) charge densities of both streams. Also the stability against radial 
perturbation is considered. Finally, in Sec. 5, we summarize our main results. 
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2 Static Relativistic Thin Disks 



In this section we present a summary of the procedure to obtain electrostatic (magne- 
tostatic) axially symmetric thin disks. The simplest metric to describe a static axially 
symmetric spacetime is the Weyl's line element 

where $ and A are functions of r and z only. The Einstein-Maxwell field equations, in 
geometrized units such that SvrG = c = ^iq = Eq = 1, are given by 

Rab = Tab, (2a) 

Tab = FacFf^'^ — -QabTcdF'^''', (2b) 
i^"';b = 0, (2C) 
Fab = Aa,b - Ab,a, (2d) 

where all symbols are understood. For the metric (|l|), the Einstein-Maxwell equations 
in presense of purely electric field are equivalent to the system 

1 

$rr + -$r + $..-^(V'' +V'') = 0, (3a) 
1p,rr + ^i^,r + i^,zz " 2($,^^/',^ + ^,z'lp,z) = 0, (3b) 



A,. = r{<l>l - $M - 'A^ii^l - i^l), (3c) 



-2* 
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A,, = 2r<l>,,$,, - re-''^^,,^,,, (3d) 
and in the magnetostatic case to 

1 

+ + - —(4 + 4) = 0, (4a) 
+ + 2(A,$,, + A,$,,) = 0, (4b) 

2$ 

A,. = r(<|.^.-$i) + ^(4-4)' (4c) 
A,, = 2r<l>,,<l>,, + ^e^^ArA., (4d) 



where ijj and A are the magnetostatic and electrostatic potencial, respectively, which 
are also functions of r and z. 

In order to obtain a solution of (0) - (0) representing a thin disc at 2; = 0, we 
assume that the components of the metric tensor are continuous across the disk, but 
their discontinuous first derivates on plane z = 0, with discontinuity functions 

Thus, the Einstein-Maxwell equations yield an energy-momentum tensor = S{z) 
and a planar current density i" = 2F"'^6{z) , where S{z) is the usual Dirac function with 
support on the disk and 

Qt = lib'^'Sb - b'^K + 9^% - 9''K + K{9''5t - 9'''5l)} 

is the distributional energy-momentum tensor. The "true" surface energy-momentum 
tensor (SEMT) of the disk, S*^, can be obtained through the relation 

S'a = j Tl ds^ = e^-* , (5) 

where dSn = ^J'gl'z dz is the "physical measure" of length in the direction normal to 

the disk, and the current density as = e^~*i'^. For the metric (|I|), the non-zero 
components of S\ are 

SI = 2e*"^{A„-2<l>„ }, (6a) 

Si = 2e*"^A„, (6b) 

and current density equal to 

= -2e*-V., (7a) 

= -2e*-^A,„ (7b) 

in the electrostatic and magnetostatic cases, respectively. All the quantities are evalua- 
ted ai z = 0"*". 

With an orthonormal tetrad ea* = {V'', W'^,X^, Y''}, where 

V = e-* (1,0,0,0) , (8a) 

W = — (0,1,0,0), (8b) 
r 

X- = e*-^(0, 0,1,0) , (8c) 

Y" = e*-^(0, 0,0,1) , (8d) 
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we can write the metric and the SEMT in the canonical forms 

9ab = -VaV, + WaW, + X,X, + Y^Y, , 



(9a) 



Sab = eVaVb+p^WaWt , (9b) 

where 

6 = -^0° , = Si , (10) 

are, respectively, the energy density and the azimuthal pressure of the disk. 



3 The Counter-Rotating Model 

We now consider, based on references |^ and [^], the possibility that the SEMT 5'°'' 
and the current density can be written as the superposition of two counter- rotating 
fluids that circulate in opposite directions; that is, we assume 

gab ^ gab ^ ^ab ^ ^^^^-^ 



(11b) 



where the quantities in the right-hand side are, respectively, the SEMT and the current 
density of the prograd and retrograd counter-rotating fluids. 

Let = {Uj_, U^, 0, 0) be the velocity vectors of the two counter-rotating fluids. In 
order to do the decomposition (|lla|) and ( p. lb|) we project the velocity vectors onto the 
tetrad ea^, using the relations |]2^ 



and thus 



With the tetrad (ISl) we can write 



ya 



^1 - U^U+f/^ - v^i - u^u_t/^ 



(12) 



(13) 



:i4a) 



1 - \]lui 



1 - \]IU'' 



(14b) 



u+ - u_ 

where U± = U]./U^ are the tangential velocities of the fluids with respect to the tetrad. 
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Using ([T^), we can write the SEMT as 



/(U_,U_)(1-U^) UlU' 
(U+-U_)2 

/(U+,U+)(1-U^J UlU^ 
(U+-U_)2 

/(u+, u_)(i - u^)5(i - \]Dh{uiut + u'^uX, 



S' 



ab 



+ 



where 



/(Ul,U2 



(u+ - u. 



eUiU2 + ■ 



(15) 



Clearly, in order to cast the SEMT in the form ( |lla|) , the mixed term must be absent 
and therefore the counter-rotating tangential velocities must be related by 



/(U+,U_ 







(16) 



where we assume that |U±| ^ 1. Then, assuming a given choice for the counter- rotating 
velocities in agreement with the above relation, we can write the SEMT as ( |lla| ) with 



so that we have two counter-rotating dust fluids with energy densities given by 

1 



(17) 



(18) 



Thus the SEMT S""^ can be written as the superposition of two counter-rotating 
dust streams if, and only if, the constraint (|T6|) admits a solution such that U+ ^ U_. 
This result is completely equivalent to the necessary and sufficient condition obtained 
in reference 



Similarly, we can write the current density in both cases as (|llb| ) with 

j± = <y±UX (19) 
where a± are the counter-rotating rest-charge densities of the fluids which are given by 



1 



Ut -u± 



U±-Ut 



(20a) 
(20b) 
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respectively. Note that the counter-rotating energy densities e± and electrostatic (mag- 
netostatic) charge densities (Te±(cTm±) are not uniquely defined by the above relations, 
also for definite values of U±. 

Another quantity related with the counter-rotating motion is the specific angular 
momentum of a particle rotating at a radius r, defined as h± = g^p^U"^. We can write 

= (21) 

This quantity can be used to analyze the stability of the disks against radial perturba- 
tions. The condition of stability, 

A_i>o, (22) 
is an extension of Rayleigh criteria of stability of a fiuid in rest in a gravitational field 
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We shall now analyze the possibility of a complete determination of the vectors f/^. 
As we can see, the constraint (|T6]) does not determine U± uniquely, and so there is 
a freedom in the choice of [/^. A possibility, commonly assumed, is to take the two 
counter-rotating fiuids as circulating along electro-geodesies 

^e±^7a6,rf/±f/i = -O^FraUl- (23) 

Let = U]_/U^ be the angular velocities of the particles. Using (|13D , ([T8|) and (|20a]) , 
in electrostatic case (^) takes the form 

5'll,r^^ + 900,r = ^^,r, (24) 

so that 

I 2 9oo,r 2j Vq 

UJ± = ± UJ , Ul = ip^r, (25) 

gu,r ^gil,r 



and similarly, using (|20b| ) instead of (|20a)) in magnetostatic case we obtain 

9ll,r^^ + 900,r = ^^,r, (26) 

so that 

_l_ 2 900,r 2j Vq , , 

U!± = ± U , UJ = A^r- (27) 

9n,r ^9ll,r 

Note that in both cases, the two geodesic fiuids circulate with equal and opposite 
velocities. 
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In order to see if the geodesic velocities agree with (p!6[), we need to compute /(U+, U_). 
In terms of uj± we get 

U± = ± U = ± 



(28) 



and so, using the Einstein-Maxwell equations (Q) and (|) and the expressions ( ]6a| ) - ( |T5D 
for the SEMT and the current density we can show that /(U+,U_) vanishes in both 
cases, so that the constraint ([16|) is equivalent to 



e 



(29) 



as is commonly assumed in the works concerning counter-rotating disks. We now have 
two counter-rotating charged dust streams with equal energy densities 



specific angular momenta 
charge densities 



h± = re 



^ Pip 



t — 



Pip' 



1 L Pv 



2r v 



(30) 
(31) 

(32a) 
(32b) 



and velocities given by (^91) . As we can see, in this case we have a complete determination 
of all the quantities involved in the CRM. 



4 Some Simple CRM Models 

4.1 CRM for Reissner-Nordstrom like disks 

The simplest electrostatic solution of the Einstein- Maxwell equations is the well-known 
Reissner-Nordstrom solution [26|, which can be written as (|l]) with 



A 



2 



In 



X + aY 
r - 1 ' 

rj.2 _ y2 



V2b 

X + a' 
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(33a) 
(33b) 
(33c) 



where a = m/k, b = e/k, with k"^ = th? — e^, so that = 1 + 6^. Here m and e are 
the mass and the charge, respectively, x and y are the prolate spheroidal coordinates, 
related to the Weyl coordinates by 

2kx = + {z + zo + kf + ^r^ + {z + zq - k^, (34a) 

2ky = y/r^ + + + ky - ^Jr'^ + (2 + - W- (34b) 

Note that we have displaced the origin of the z axis in zq. This solution can be genera- 
ted, in these coordinates, using the well-known complex potencial formalism proposed 



by Ernst [p7| , |28[| from Schwarzschild solution p9[. Indeed for 6 = 0, it goes over into 
the Schwarzschild solution. Thus, h is the parameter governing the electric field. One 
can obtain its magnetostatic counterpart cumputing the magnetostatic potencial via 

A, =kf-\\-y')i,^y, (35a) 

Ay =-kf-\x' (35b) 

Thus, we find 

A = V2kby, (36) 

again with 0? — iP' = 1. From the above expressions we can compute the physical 
quantities associated with disk. We obtain 

. ^ 4y{x^ - l){ax + f) 

^ (x + a)2(x2 -y2)3/2' I ) 

4xy{l — y^) 

= (a; + a)(x2-y2)3/2' (38) 

L = (39) 

(x^ — y'^y/'^{x + a)^' 



2V2bx(l 



^'^ (x2-y2)\/2(/+a)' (40) 

where e = ke, = kp^ and jt = /cj*. x and y are given by 

2x = v'r^ + (a + 1)^ + y^r^ + (a - 1)2, (41a) 



2y = + {a + 1)2 - v/r2 + (a - 1)2, (41b) 
where r = r/k and a = zo/k, with a > 1. 
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In order to study the behavior of these quantities we perform a graphical anahsis of 
them for disks with a = 1.4 and 6 = 0, 0.5, 1.0, and 1.5. In Fig. |l](a) we show the 
energy density e as a function of f. We see that the energy density present a maximum 
at f = and then decreases rapidly with f. We also see that the presence of electric 
(magnetic) field decreases the energy density in the central regions of the disks and later 
increases it. For other values of a and b we obtain a similar behavior. 

In Fig. |l](6) we plot the azimuthal pressure as function of f. We can observe 
that the pressure increases rapidly as one moves away from the disk center, reaches a 
maximum and later rapidly decreases. Note that electric (magnetic) field decreases the 
pressure everywhere on the disk. For other values of a and b we have a similar behavior. 

The electrostatic current is presented in Fig. @(a). Similarly to the energy density, 
it has a maximum at the disk center and then decreases rapidly with f. For other values 
of a and b we observe a similar behavior. In addition, the magnetostatic current j;^ is 
ploted in Fig. ^6) as function also of f. As we can observe it has a similar behavior to 
the pressure. 

We now consider the CRM for the same value of the parameters. All the significant 
quantities can also be expresed in analytic form from above expressions but the results 
are so cumbersome that it is best just to analyze them graphically. In Fig. ^(a) we 
show the velocity curves of counter-rotating streams as functions of r. We observe 
that it increases rapidly in the central region of the disk, achieves a maximum and later 
decreases monotonly. We also see that the presence of electric (magnetic) field makes 
least relativistic the disks. In addition, in Fig. ^(6), we plot for disks with b = 0.5 
and a = 1.01, 1.1, 1.4, and 2.0. We find that the disks become least relativistic with 
increasing a. We also find that the disks with a < 1 cannot be built from CRM because 

> 1 (not shown in the figure). 

In Fig. ^(a) we have drawn the specific angular momentum h"^ of counter-rotating 
fluids, where h = h±/k. In the cases considered we obtain h"^ as increasing monotonic 
functions of f what correspond to stable CRM for the disks. However, the CRM cannot 
be applied for 6 = 4.0 (dotted curve). Thus the presence of electric (magnetic) field 
can makes unstable the CRM against radial perturbations. Finally, in Figs. ^(6) and 
H the plots of the mass densities e± and electrostatic (magnetostatic) charge densities 
c"e±(c"m±) of both streams are shown. These present a maximun at the disks center 
and then decrease monotonly. Therefore, the CRM constructed from this value of the 
parameters are well behaved. 
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4.2 CRM for Taub-NUT like disks 



A Taub-NUT like solution to the Einstein- Maxwell equations is 



A 



In 



x"^ + 2ax + 1 



2 In 



X 



x^ 



2^2 6; 



•X 



+ 2ax + 1 ■ 



(42a) 
(42b) 
(42c) 



with — 6^ = 1. X and y are, again, the prolate spheroidal coordinates, given by Eqs. 
( P4a|) and ( P4b|) . This solution can also be generated, in these coordinates, using the 
well-known complex potencial formalism proposed by Ernst |2S| from Weyl 2-solution 
(Darmois ||26|). Indeed for 6 = 0, it goes over into the Darmois solution . Thus, again, b 
is the electric parameter. One can also obtain its magnetostatic counterpart cumputing 
the magnetostatic potencial using (|35|), and we get 



A = 2V2kby, 



(43) 



again with — b"^ = 1. The physical quantities associated with the disk now can be 
written as 

8y[x{ax^ — Sax — 2) + y'^{2x^ + 3ax^ — a)] 



it 



(x2-l)(x2 + 2ax + l)2 

_ IQxyjl - f) 

" (x2-l)(x2 + 2ax + l)' 

4^/2by{x^ - l){x^ - y^) 
(x2 + 2ax + 1)3 ' 

AV2bx{l - y^){x^ - f) 



(44) 
(45) 
(46) 
(47) 



(x2 - l)(x2 + 2ax + 1) ' 

where x and y are given by Eqs. ( |41a|) and ( [41b| ). 

In Figs. H and |^ the plots of the quantities e, p^, and are presented for disks with 
a = 2.5 and 6 = 0, 0.5, 1.0, and 1.5. We see that these disks have a similar behavior to 
the previous case. Equally, the relevant quantities of the CRM are shown in following 
figures for the same value of the parameters. Theses CRM are more relativistic than 
the ones built from Reissner-Nordstrom like solution (Fig. ^(a) ). Note that (Fig. ^(b) 
) the disks with b = 0.5 a = 1.5 (solid curve) cannot be constructed from the CRM 
because > 1. We also note that the presence of electric (magnetic) field can makes 
unstable the CRM against radial perturbations (Fig. ^(a)). Thus the CRM cannot 
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apply for 6 = 3 (dotted curve) . The remaining functions are plotted for the value of the 
parameters representing a physically acceptable CRM (Fig. ^(6) and |l^). Also these 
have a similar behavior to the previous case. 



4.3 CRM for Kerr like disks 

A Kerr-like solution to the Einstein-Maxwell equations is 



A 



In 



2 In 



[ax 



1)2 _ 52^2 

- hW - ] 



a?{x'^ — 2/2) 
2^/2hy 



{ax + 1)2 — 62y2 • 



(48a) 
(48b) 
(48c) 



with a^ — IP' = 1. X and y are, again, the prolate spheroidal coordinates, given by Eqs. 
( ^4a| ) and ( p4bD . This solution can be generated, in these coordinates, using a well-know 



theorem proposed by Bonnor I^O], ^ from Kerr solution. The previous solution can also 
be obtained using this same theorem from Taub-NUT solution. For 6 = 0, it also goes 
over into the Darmois solution. Thus, again, h is the electric parameter. One can also 
obtain its magnetostatic counterpart cumputing the magnetostatic potencial using (p5D, 
and gives 

V2kh{l-y^){ax + l) ^ 
ai^a^x"^ — h'^y'^ — 1) ' 

again with a^ — IP' = 1. This solution describes the field of a massive magnetic dipole. 
The physical quantities associated with disk now can be written as 

e = 8a4y{(x2-^2)[^^^2_^)[^^^^;L)2 + 6V]-262a;(aa; + l)(l-y2)] 

-2x{x^ - 1)(1 - y^)[{ax + 1)^ - feV]} (50) 
/{a^x^ - b^f - l)2[(ax + 1)2 - b^f, 



lQa'^xy{x'^ — 1)(1 — y"^ 



(a^x^ — 52^2 _ \Y\J^ax + 1)2 — b'^y 



i721 ' 



(51) 



4V2a%{x^ - y^){x{l - y^)[{ax + if + ^2^2] _ 2ay\ax + l){x^ - 1)} 



(a'^x'^ 



b^f-l)[{ax + iy-b^y 



,-7213 



2V2a%y{x^ - 1)(1 - - y^)[iax + l)(3ax + 1) + b^f] 

(a2x2 — 52^2 _ 1)3|^('q^2; -|- 1)2 _ 52|^2j 



(52) 
(53) 
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where x and y are given by Eqs. ( |41a| ) and ( [41b| ). 

In figs. |Tl| and the plots of the physical quantities describing the disks are shown 
for a = 1.8 and 6 = 0, 0.5, 1.0, and 1.5. The energy density behaves in the opposite 
way to the previous cases. That is, near to the disk center it increases when the electric 
(magnetic) field is applied and decreases later. The other quantities have a similar 
behavior to the precedent cases. However, the electrostatic current after some r takes 
negative values. Likewise, the quantities corresponding to the CRM are shown in follo- 
wing figures for the same value of the of the parameters a and b. Theses CRM are more 
relativistic than the ones above considered (Fig. |l3|(a)). Note that (Fig. |l3|(b)) the disks 
with b = 0.5 and a = 1.4 (solid curve) cannot be constructed from the CRM because 

> 1. We also note that the presence of electric (magnetic) field can stabilize the 
CRM against radial perturbations (Fig. |I^(a)). Therefore, only the CRM constructed 
with b = 1.0, and 1.5 are well behaved. Finally, the other functions are drawn for the 
same value of the parameters (Figs. |T^(&) and[T^). 



5 Discussion 

A detailed study of the Counter-Rotating Model for generic electrostatic (magnetos- 
tatic) axially symmetric thin disks without radial pressure was presented. A general 
constraint over the counter-rotating tangential velocities was found, needed to cast the 
surface energy-momentum tensor of the disk in such a way that it can be interpreted 
as the superposition of two counter-rotating charged dust fluids. The constraint found 
is completely equivalent to the necessary and sufficient condition obtained in reference 
p3| . We next showed that this constraint is satisfied if we take the two counter- rotating 
fiuids as circulating along electrogeodesics with equal and opposite tangential velocities. 
We also have obtained explicit expressions for the energy densities, electrostatic (mag- 
netostatic) current densities and velocities of the counter-rotating streams in terms of 
the energy density, azimuthal pressure and planar current density of the disk, that are 
also equivalent to the correspondig expressions in reference p3| . 

Three specific examples were considered in the present work based in simple solu- 
tions to the Einsteins-Maxwell equations generated by conventional solution-generating 
tecniques [26|. We found that the CRM for Kerr-like disks are more relativistics than 



the ones obtained from Taub-NUT and Reissner-Nordstrom like solutions. We also saw 
that the presence of electric (magnetic) field can make unstable the CRM against radial 
perturbations in the case of Taub-NUT and Reissner-Nordstrom like disks, and con- 
versely, stabilize the CRM in the case of Kerr-like disks. We also constructed some 
CRM with well defined counter-rotating tangential velocities and stable against radial 
perturbations. 

On the other hand, the generalization of the Counter- Rotating Model presented here 
to the case with radial pressure is in consideration. Also, the generalization to rotating 
thin disks with or without radial pressure in presence of electromagnetic fields is being 
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considered. 
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Figure 1: (a) Energy density e and (6) azimuthal pressure as functions of f for disks 
with a = 1.5 and 6 = (solid curve), 0.5, 1.0, and 1.5 (dotted curve). 



r r 
(a) (6) 

Figure 2: Surface current density: (a) and (6) as functions of f for disks with 
Q! = 1.5 and 6 = (axis f), 0.5, 1.0, and 1.5 (dotted curve). 
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Figure 3: Tangential velocity as function of r for disks with (a) a = 1.5 and 6 = 
(solid curve), 0.5, 1.0, and 1.5 (dotted curve), and (b) b — 0.5 and a — 1.01 (solid curve), 
1.1, 1.4, 2.0 (dotted curve) 




Figure 4: (a) Specific angular momentum /i^ as function of f for disks with a = 1.5 and 
b — (sohd curve), 0.5, 1.0, 1.5 and 4.0 (dotted curve), (b) Mass densities e± as function 
of f for disks with a — 1.5 and b — (solid curve), 0.5, 1.0, and 1.5 (dotted curve). 
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Figure 5: Charge densities: (a) ae± and (6) am± as functions of f for disks with a — 1.5 
and 6 = (axis f), 0.5, 1.0, and 1.5 (dotted curve). 




Figure 6: (a) Energy density e and (6) azimuthal pressure as functions of f for disks 
with a = 2.5 and 6 = (sohd curve), 0.5, 1.0, and 1.5 (dotted curve). 



18 



r r 
(a) (b) 

Figure 7: Planar current density: (a) it and (6) ji^ as functions of f for disks with a — 2.5 
and 6 = (axis f), 0.5, 1.0, and 1.5 (dotted curve). 




Figure 8: Tangential velocity as function of f for disks with (a) a — 2.5 and 6 = 
(solid curve), 0.5, 1.0, and 1.5 (dotted curve), and (6) b — 0.5 and a — 1.5 (solid curve), 
1.7, 2.0, 2.5 (dotted curve) 
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(a) (b) 

Figure 9: (a) Specific angular momentum h"^ as function of f for disks with a = 2.5 
and 6 = (solid curve), 0.5, 1.0, 1.5, and 3.0 (dotted curve). (6) Mass densities e± as 
function of f for disks with a — 2.5 and b — (solid curve), 0.5, 1.0, and 1.5 (dotted 
curve) . 
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Figure 10: Charge densities: (a) ae± and (6) cr^-t as functions of f for disks with a — 2.5 
and b — (axis f), 0.5, 1.0, and 1.5 (dotted curve). 
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Figure 12: Current density: (a) and (6) }^ as functions of f for disks with a — 1.8 and 
6 = (axis f), 0.5, 1.0, and 1.5 (dotted curve). 
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Figure 13: Tangential velocity as function of f for disks with (a) a = 1.8 and 6 = 
(solid curve), 0.5, 1.0, and 1.5 (dotted curve), and (6) b = 0.5 and a = 1.4 (solid curve), 
1.6, 2.0, and 2.5 (dotted curve) 




Figure 14: (a) Specific angular momentum h'^ and (6) mass densities e± as functions of 
f for disks with a — 1.8 and b — (solid curve), 0.5, 1.0, and 1.5 (dotted curve). 
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Figure 15: Charge densities: (a) ae± and (6)crTO± as functions of f for disks with a — 1.8 
and b — (axis f), 0.5, 1.0, and 1.5 (dotted curve). 
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